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We study the effects of braiding interactions on single anyon dynamics using a quantum walk 
model on a quasi- 1-dimensional ladder filled with stationary anyons. The model includes loss of 
information of the coin and nonlocal fusion degrees of freedom on every second time step, such that 
the entanglement between the position states and the exponentially growing auxiliary degrees of 
freedom is lost. The computational complexity of numerical calculations reduces drastically from 
the fully coherent anyonic quantum walk model, allowing for relatively long simulations for anyons 
which are spin— 1/2 irreps of SU(2)fc Chern-Simons theory. We find that for Abelian anyons, the 
walk retains the ballistic spreading velocity just like particles with trivial braiding statistics. For 
non-Abelian anyons, the numerical results indicate that the spreading velocity is linearly dependent 
on the number of time steps. By approximating the Kraus generators of the time evolution map 
by circulant matrices, it is shown that the spatial probability distribution for the k = 2 walk, 
corresponding to Ising model anyons, is equal to the classical unbiased random walk distribution. 



I. INTRODUCTION 



Quantum walks have proven to be a fruitful model to study the dynamics of particles on lattices or more general 



graphs. The effects of exchange statistics on multiparticle quantum walks has been studied for bosons and fcrmions 
with a notable effect of statistics. For example it has been shown that quantum walks with non-interacting bosons 
give rise to effective interactions which lead to Bose-Einstein condensation in graphs with spatial dimension d < 2 
[1]. Quantum walks with pairs of bosons or fermions have been studied in the context of the graph isomorphism 
problem by relating the Green's function of the evolution to the spectrum of the graph. It was shown that these 

,_^ models have much more computational power when the particles are allowed to interact because they can distinguish 

^ non- isomorphic graphs that noninteracting particles cannot [ ]. Recently using integrated photonics there has been 

OS experimental realization of two particle quantum walks which simulated bosonic, fermionic, and semionic (a type of 

OS Abelian anyon) statistics [ ]. 

A natural extension of multiparticle quantum walks is to consider anyons which are particles that exist in two 
dimensions and have richer exchange statistics than bosons or fermions. In the case of Abelian anyons, exchanges, or 
braids, yield some model dependent phase e 1 ^ while for non-Abelian anyons there is a large Hilbert space of fusion 
degrees of freedom of the multi particle system and a braid performs a unitary transformation on this vector space. 
A quasi-1-dimensional model for anyonic quantum walks was introduced in Ref. [ ] where it was shown that Abelian 
anyons have quadratic variance just like standard quantum walks with trivial statistics but that non-Abelian anyons 
create entanglement during the walk purely due to braiding that acts to decohere the spatial degree of freedom. In 
Ref. [.}] the quantum walk on a quasi- 1-dimensional ladder was studied in detail for one type of non-Abelian anyon 
known as the Ising model anyon. This anyon corresponds to a spin— 1/2 irrep. of SU(2)k=2 Chern-Simons theory, 
and the variance of the probability distribution of the walker was proven to be asymptotically linear in time with 
coefficient one, corresponding to the classical random walk behaviour. That particular case points out that there is a 
fundamental difference between the dynamics of Abelian vs. non-Abelian anyons. 

In this work we explore these differences for more general anyon models, specifically for spin— 1/2 irreps of all 
values of the index k in SU(2)k Chern-Simons theory. Numerical results show that for small number of time steps, 
the position distribution approaches the quantum walk distribution as index k grows [ ], and it was conjectured that 
the distribution looks classical when fc <C t. The problem with the fully coherent anyonic quantum walk is that to 
evaluate the probability distribution for t iterations of the discrete evolution operator, the number of paths which 
contribute to the walk grows exponentially with t. It turns out that evaluating each path is equivalent to computing 
the Jones polynomial for the link corresponding to the closed world line of that path, a computation which is itself 
exponentially hard (except for the cases k — 1,2,4 [(>]) in the number of time steps t. To simplify the calculations 
for general k, we adopt here a special quantum walk protocol which introduces loss of memory to the coin and fusion 
degrees of freedom. Such a protocol allows for calculating the walk distribution for a relatively large number of time 
steps and to use certain approximations to diagonalize the walk operator and obtain analytical results. 



II. ANYONIC V 2 QUANTUM WALK 

The anyonic quantum walk is a dynamical model for a single anyon (mobile or "walker" anyon) moving in a lattice 
of fixed anyons (stationary "background" anyons) . The mobile anyon has an extra degree of freedom called the coin 
which is coupled to the direction of its movement. Here we consider a quasi- 1-dimensional lattice of background anyons 
where the walker has only two directions of movement, left or right, and the coin is thus two-dimensional. Although 
the anyons are arranged on a line, they exist on a two-dimensional manifold such that the mobile anyon can pass 
them from above or below without coming into contact with the stationary anyons. This dynamics is accomodated 
by a walk on a ladder where the top(bottom) leg corresponds to mode |0)(|l)) of the coin, the sites correspond to 
the rungs and are labelled by an integer s, the stationary anyons are pinned in the islands of the ladder (see Fig. 1). 
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FIG. 1: The anyonic quantum walk on a ladder of point contacts. The red oval is the walker anyon which moves along the 
black lines, and the grey ovals are the stationary anyons which occupy the islands between the point contacts. The anyons are 
initialized in the state |$o) meaning they are created out of the vacuum in pairs as indicated by the ovals joined by strings, 
and half of each pair is outside the ladder not participating in the walk. When on the upper(lower) edge the walker always 
to the left (right). The occupation sites of the mobile anyon are labelled by the positions s of the point contacts where the 
mobile anyon can tunnel. The tunneling matrix between the legs at each rung corresponds to the coin shuffling operator Uy- 
Each lattice shift corresponds to a braid between the worldlines of the walker anyon and a stationary anyon, represented by 
the generator b s . 



The arrangement in Fig. 1 is a crude model of quasiparticle dynamics in Fractional Quantum Hall samples which 
support anyonic quasiparticle excitations [7, 8]. The edge modes (horizontal in Fig. 1) carry mobile quasiparticles 
with a current direction determined by an externally applied voltage, and bulk quasiparticles can be pinned at antidots 
within the islands. The upper and lower edge modes are mixed by introducing point contacts controlled by front gates 
between the islands, allowing tunneling from edge to edge. The coin is thus encoded as propagation to left or right 
on the upper or lower edges respectively, and the scattering matrix that describes the tunneling between the edges 
represents the quantum coin flip operator. The model assumes a hardcore approximation for the anyons, such that 
they are always far apart, no polarization to fusion channels occurs and the anyons interact only non- locally via their 
braiding statistics. The model is also chiral: the particles on the upper edge are only allowed to move left and on the 
lower edge only to the right, and all braids are therefore anti-clockwise. 

When the walker shifts between lattice sites, the wave function is multiplied by the braid generator that corresponds 
to braiding the walker and the stationary anyon. The total Hilbert space is the tensor product space of the spatial 

sites, fusion degrees of freedom and the coin, % — T-L s -pa,cc®T~i- fusion •Sf'Hcoin- The position space is spanned by { \s) } s=0 , 
the fusion space accommodates 2N anyons (N — 1 stationary anyons, the walker anyon and their antiparticle pairs) 
with total charge vacuum, and the coin space is spanned by { 0), l)}- We assume periodic boundary conditions, and 
choose the number of sites to be always larger than the length of the walk, N > At + 1. The walker starts initially 
from the localized state So)(so|- For reasons described below, rather than fully coherently evolving the system over 
many steps, the walk evolves coherently for two cycles followed by tracing over fusion and coin degrees of freedom. 
Afterward the fusion state is reset to the vacuum pair state $o/(^0i the coin state to Co)(co , and the dynamics is 
iterated: 



Ps (t + 1) = Tr fjC ((£/ s [/ F ) 2 Ps (t) ® |$ )($ | <g> |c >(co| ((IWf) 2 ) 1 ) 
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U s = J2 { S s+i ®b s ®P + Sf ® b s ® Pi) (4) 

when the Hadamard matrix _ff is chosen as the coin flip operator, Sj" = s — l)(s and 5+ = s + l)(s are the shift 
operators and Po = |0)(0|, Pi = l/(l are projectors to the coin states. The braid generators b s are the matrices 
associated with braiding the anyons s and s + 1 in anticlockwise direction. Due to the periodic boundary conditions, 
there is also a braid generator &jv_i associated with braiding the anyons 1 and N which can be expressed in terms of 
the other generators 1 , but we choose the system size so that the walker never crosses the boundary and this generator 
is never applied. A string of braid generators, B = Y[ s . b s . is called a braid word. The braid matrices are determined 
by the anyon model in question. In a physical setting, the properties of the underlying material supporting the anyonic 
excitations determine the anyon model. 

We consider a general class of non-Abelian anyons corresponding to spin— 1/2 irreps of the quantum group SU(2)fc 
which is a group with a deformed version of the SU(2) algebra. The fusion rules of the SU(2)fc model satisfy the 
triangle inequality and integer sum condition as in SU(2): 

jl+32 

h x h = Yl J 

3=\h-h\ 
but with two restrictions on the total spin charge j: 

3 < k/2] h +32+3 <k 

where j\ and j-z are two individual charges and the parameter k is the level of the theory. The Hilbert space of SU(2)fc 
anyons grows as dim(Hf us ; on ) ex d N where d = 2cos{-j^) is the quantum dimension. Since the fusion space grows 
exponentially in the total number of anyons, if one is keeping track of correlations of the walker with the other anyons, 
the number of degrees of freedom to keep track of grows very fast with the number of time steps and the calculations 
become inefficient. To calculate the probability distribution of the walker for a given time step, only the calculation of 
the trace over the fusion degrees of freedom is necessary. Fortunately, there exists a connection between the trace over 
the representations of the braid group <$o \B ^o) and link invariants from knot theory that simplifies calculations of 
the trace. The trace over the fusion degrees of freedom can be understood in two equivalent ways: the usual matrix 
trace over the bracket representation of the density matrix of the fusion state, or the diagrammatic quantum trace 
which is defined by taking the fusion diagram representation of the state and joining the open ends of the incoming 
and outgoing charge indices [ ]. The choice of which indices are joined together determines a tracing scheme for the 
braid word, and the tracing corresponds to a closure of a braid in knot theory, such that the closed braid forms a link. 
The usual tracing schemes are plat tracing, where the neighbouring charges are fused such that their total charge is 
vacuum, and Markov tracing, where each charge is connected to itself before and after time evolution. 

We choose the initial state of the anyons |$o)($o| such that the tracing corresponds to Markov tracing. The anyons 
are initially created from vacuum pairs, after which one member of each pair is dragged out of the system. In a braid 
presentation, this corresponds to having nearest-neighbour pairs with total charge vacuum, followed by a braid word 
which moves all left members of the pairs to the left. The total number of background anyons in an iV-site walk is 
N—l, but the total Hilbert space accommodates 2N anyons (background anyons, walker anyon, and their antiparticle 
pairs) . Only half of the anyons participate in the quantum walk and after the two-step evolution all anyons are fused 
back together with their antiparticle pairs, see Fig. 2. In this scheme, the expectation value ($o|-5 n&o) corresponds 
diagrammatically to the Markov trace over the braid word B. 

The virtue of using the state |$o) as the initial fusion state is that the expectation value ($o|-^|^o) can now be 
expressed in terms of the link invariant called the Kauffman bracket which is a Laurent polynomial in the variable A 
[10], denoted (L}(A), where L is a link, or a collection of tangled loops. This is written formally as 

<*o|B|*o> = ^ (5) 

where L = (5 Markov ) [ s the link obtained from the braid word B via closure by the initial state \&o), d is the 
quantum dimension of the anyons, n is the number of distinct components (strands) in the link L, and the value of 
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FIG. 2: A link generated by the anyon worldlines in a quantum walk. Particle-antiparticle pairs are initially created from 
the vacuum state and one member of each pair is dragged out of the system. The walker anyon braids with the background 
anyons for two steps, after which each particle is fused back with their antiparticle pairs. Such a protocol implies that the 
diagrammatic tracing scheme is Markov. 

the parameter for SU(2)fc anyons is A = ie~ l7T / 2 ( k+2 \ The calculation of the trace over the fusion states thus reduces 
to evaluation of the Kauffman bracket of the link that is drawn when the walker anyon moves on a particular path. 
The evaluation of the Kauffman bracket is not necessarily more efficient than calculating the matrix trace of the braid 
representations, but if one restricts to a model where the walker evolves coherently only for a finite amount of time 
steps, it turns out that for arbitrary long walks only a fixed number of links needs to be evaluated and the quantum 
walk distribution can be calculated exactly for fairly large number of time steps. 

In the usual quantum walk protocol, the quantum speedup occurs because the position of the walker becomes 
entangled with the coin. The system evolves coherently and the correlations between the position states and coin 
states preserve the memory in the system, hence the system dynamics is highly non-Markovian. If the walk is subject 
to decoherence, the correlations between the position and the coin become degraded and some of the memory in the 
system is lost. The loss of memory can be modeled using the so called u V n model" [ ], where the state of the coin is 
erased and reset on every nth step, for example by doing a measurement on the coin and preparing it over again in the 
initial state. After resetting the coin, the system is in a product state and the correlations between the position and 
coin are lost. Remarkably even in the V 2 model there is enough coherence left in the system to provide for quadratic 
speed up over the classical random walk. Specifically, it was shown in Ref. [ '.] that the variance asymptotically scales 
like cr(t) 2 — K2t 2 + K^t for some constants K2,K^. 

Here the dynamics of our anyonic walk, Eq. 1, is the V 2 model where V = UsUf- We write the dynamics as a 
super operator on the density operator p s (t) for the spatial degrees of freedom of the walker: 

Ps (t + 1) = S(p s (t)) 



E (l®(f\®(c\)(U S U F ) 2 A,(t)®|#o)(*o|®|cto)(co| ((U s U F ) 2 y(l®\f)®\c) 



f,c 



E {l®{f\®{c\)(U s U ¥ ) (r®|* )®|co» p B (t) 

f,c L J 



{l®(<S> \®(c \)({U s U F ) 2 y(l®\f)®\c)) 
= EE f cPs(t)E} c 

f,c 

where the Kraus generators have been defined as 

E fc = (I® (f\ <8> <c|) {U S U F ) 2 (J® |* > ® |co»- 
Using equations (1)— (4) we find for the double step operator 

(U S U F ) 2 = E [\s-2)(s\<S>b s - 2 b s - 1 ®PoHPoH + \s)(s\ ® b 2 ® P HP X H 



(6) 



(7) 



+ \s)(s\ <g> 6^_ 1 <g> P X HP H + \s + 2)(s| ® b s+1 b s ® PxHPxH 



It is clear from this equation that the double step operator shifts the coefficients of the density matrix by two 
rows or columns, or not at all. By inserting this expression to the superoperator and using the completeness 

of the fusion and coin bases, £ (/|6»|*o)(*o|&]k|/> = (*o|&J/&»|*o) and Y,(c\PaHP b H\co)(co\WP a >H^P b ,\c) = 

/ 
3a,b' (co\H^ P a > H^ P a HPbH |co), the action of the superoperator on a general element |s)(s'| of the spatial density 
matrix can be written as a sum of 7 terms: 

£ (\ s )( s '\) = 3[|s-2)( S '-2|($ |6^ 1 6 t s ,_ 2 6 s _ 2 6 s _ 1 |$o) + \s-2)(s'\($o\btX-2b s - 1 \$ ) 

+ | s )< s '-2|($ |& t s ,_ 1 6 t s ,_ 2 62|$ )+ | s )( 5 '| ((* |&J?&2|*o) + (*o|&I?-i&2-i|*o>) 

(8) 

- \s)(s' + 2\ (*o|&J,&J/ +1 6 2-i|*o) - |s + 2>(s'| (* |&t?-i & «+i&«l*o) 

- | s + 2>(s' + 2| ($ \bl,bl, +l b s+1 b s \$ ) 

where we have chosen \cq) = |0) such that (co\H^ P a i H^ P a H P b H\co) = ±|. 

Equation (8) shows that for any element of the initial spatial density matrix, there are only seven nonzero elements 
of the superoperator £ . Writing the density matrix in a vectorized form, the superoperator can be written as a matrix: 
p' = ^2 f c (Efc <£> Ef c )p> wnere the superoperator matrix Y^f c (Efc ® Ef c ) is a band matrix with seven diagonal 
bands, and the rest of the elements are zero. Starting from the initial position density matrix p s (0) — So)(so , the 
position density matrix after t time steps is obtained by applying the superoperator t times: 

A,(*)=f*(ft(0)). (9) 

The elements of the superoperator are given by the product of the coin term (^cq\H^ P a i H^ P a H P b H |c ) = ±| and 
the fusion term (& \B(s, s')|$o)- The probability distribution at time step t is obtained by starting from an initially 
localized state So)( s o and applying the superoperator t times. The site probabilities are then given by the diagonal 
elements of the spatial density matrix. We have calculated the variance, a(t) 2 = (s 2 ) — (s) 2 , of the V 2 walk for various 
choices of the level k up to t — 100 iterations of the superoperator, corresponding to 200 quantum walk steps. 

One interesting feature of the V 2 anyonic walk is that the effect of braiding statistics is trivial for Abelian anyons. 
If all the stationary anyons are of the same type, the walker always picks the same phase e 1 ^' 2 when braiding with 
them, such that 

(*o|B|*o) = e-^/V^V^V^ 2 = 1. (10) 

This is similar to the fully coherent anyonic quantum walk, where the wave function always picks the phase e 1 ^ 1 ! 2 
during the forward time evolution and phase e -2 **/ 2 during backward time evolution, such that the overall effect is 
trivial. Another way to think about the Abelian walk is to look at the case k = 1 in Chern-Simons theory. There are 
only two charges {0, |} in this model, the fusion channels are unambiguous: j-± x j 2 = j, and the quantum dimension 
is d = 1, so this model is Abelian. Substituting A — ie^™/ 6 to the values of the brackets in Table I in Appendix A 
gives ($o|-B( s j s ')|^o) = 1 f° r all brackets, confirming that the Abelian walk gives the same dynamics as the original 
V 2 quantum walk. The variance of the V 2 walk for Abelian anyons is plotted in Fig. 3 a). It shows the known fact [11] 
that the original V 2 walk propagates slower than the fully coherent quantum walk, but qualitatively the behaviour 
is still the same: the variance depends quadratically on time. The best fit for the variance of the k — 1 walk is given 
by a 2 =0.125 1 2 + 0.75 1. 

In the non- Abelian case, the construction of the superoperator £ requires evaluation of all eight expectation values 
in Eq. (8) for all s and s' . The evaluation can be done by using the relation between the expectation values and the 
Kauffman bracket in Eq. (5). Each braid word consists of four braid generators, so that each planar representation 
of the corresponding link has four crossings. A very convenient feature of the link diagrams is that if s and s' are far 
apart, the forward and backward braid words act on disjoint sets of strands, and the Kauffman brackets of all disjoint 
braids are equal (see Fig. 4 in Appendix A). Thus, the calculation of each expectation value involves only calculating 
the disjoint case and a small number of cases where s' is close to s. The values of the expectation values are tabulated 
in Appendix A for all unique choices of s' — s. 

The variances for various choices of k are plotted in Fig. 3 b). We observe that for k = 2, the variance follows 
exactly the linear line a 2 — t. It is interesting to note that this behaviour is exactly the same as was shown for the 
fully coherent walk with k = 2 [ >]. For the rest of the values of k, there is surprisingly small difference in the variance, 
with the slope of the variance slightly increasing as a function of k. For k = 3 and k = 4 (not shown in the figure) we 
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FIG. 3: The variance a of the two-step walk as a function of time t for various anyon models indexed by k. The scaling of 
the sites and time steps is chosen such that one iteration of the superoperator £ corresponds to two iterations of the quantum 
walk time evolution operator, and at each iteration of £ the walker shifts twice, a) Case k — 1 corresponds to the Abelian 
walk, which also coincides with the original V 2 quantum walk. Plotted for comparison are the variance of the classical random 
walk (RW) and the ordinary quantum walk (QW) with the same initial coin state and coin flip matrix, b) Non-Abelian SU(2)fe 
models zoomed to the last 100 time steps. Also plotted is the linear curve a 2 — t which corresponds to the classical random 
walk. 

find that the slope of the variance is slightly smaller than 1, k = 3 having the smallest slope (best fit 0.9877). The 
highest variance is obtained when k — > oo, which corresponds to choice of the parameters A = i and d = 2, with slope 
1.0665. It should be noted that in the fully coherent anyonic quantum walk model, taking the limit k — > oo leads 
to the conventional SU(2) algebra, and the braiding generators are a non- Abelian representation of the permutation 
group. However when we introduced decoherence in the system, the braid words that contribute to the walk have 
different structure than in the fully coherent model. In the fully coherent model, the diagonal elements of the spatial 
degree of freedom of the walker correspond to a sum over all paths where the walker strand and necessarily all other 
strands return their initial position after forward and backward time evolution. Hence each of these paths is a trivial 
permutation in the fusion space and has no effect on the diagonal elements of the quantum walk. In the V 2 model 
however, the links do not have a dedicated walker strand as can be seen from Fig. 4 in Appendix A. The nontrivial 
character of the k — > oo model thus follows from the fact that the links are not identity permutations and hence 
expectation values in Table I are not equal to 1, as they would be in the fully coherent walk. The key difference is 
that only for a highly non-Markovian environment does the k — ¥ oo case reduce to the normal quantum walk for only 
if one includes memory of all previous steps does the fusion DOF become disentangled with the spatial DOF. 



III. APPROXIMATION BY CIRCULANT MATRICES 



In the previous section, we derived the exact superoperator that describes how the spatial density matrix trans- 
forms in a single step of the walk, and obtained numerical results by applying the superoperator to the density matrix 
repeatedly. We did not however obtain any expression for the density matrix after evolution by an arbitrary number 
of steps. In the following, we approximate the Kraus operators Ef c by circulant matrices which have uniform coef- 
ficients on each diagonal, and diagonalize the Kraus operators via Fourier transform to find arbitrary powers of the 
superoperator £ in a compact form. 

Calling Eqs. (6) and (7) from previous section, the Kraus generators are written explicitly as 



Efc = £ [ C °c° (f\bsb s+ i\^o) \s)(s + 2\ + (/| {Cfbl + Cl%U) |* > |*><* 

S 

+ C c n (/|6 s+1 6 s |$ ) \s + 2)(i 



(11) 



where the coin term has been defined as C° b := (c|P a i?Pf,i7|co). 



The matrix expression 


of the generators is 






f d/c(0) 


d f c(l) 


a/c(0) 


E fc = 


bfc(0) 

a fc (N - 2) 


6/c(l) 


d /c (2) 



0/c(l) 



& /c (iV-2) 



o/c(JV - 1) 



6/ c (JV-3) 



6 /c (JV - 1) 

a /c (iV-3) 
d fc (N-l)J 



(12) 



where 



o /o («) = c c 00 </|& s 6 s+1 |a> ), 

bfc(s) = <7 C U (/|& S+1 & S |$ >, 

d/c(«) = a c 01 </|6^o) + C c 10 (/|6^ 1 |cf> ). 



(13) 

(14) 
(15) 



Comparing the matrix above with the general circulant matrix given in Eq. (Bl) of the Appendix B, we see 
that although matrix Ef c is not a circulant matrix, it can be turned to a circulant matrix if the s-dependence of 
parameters af c (s), bf c (s), df c (s) is suppressed in some meaningful way, so that the sequences of these parametrers 

get approximated by some respective sequences a/ c , bf c and df c . This suppression of the s-dependence also leads to 
circulant form for Ef c . The advantage of having the Kraus generators being circulant matrices stems from the fact 
that the spectral decomposition problem of such matrices is solved (see Appendix B), and therefore this enables us 
to express the CP map of our QW and its powers in terms of the orthogonal eigenprojections of the generators. To 
this end we introduce for each a/ c (s), bf c (s) and df c (s) the following approximation for each s G 7Ln 



a fc (s) 
dfc(s) 
b fc(s) 



a fc (s) 
dfc(s) 



1 

N 
1 

N 
1 



Tr(h 2 E fc ) 
Tr(E fc ) 



bfe(s) = jjTr(h- 2 E fc ) 



where the elementary band matrices are defined as h = ^2 se 

N ^ 
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(f\b s b s+1 \$ ), 
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(16) 
(17) 
(18) 

(19) 
(20) 
(21) 



i.e. the parameters a,f c (s), bf c (s), df c (s) are approximated by the arithmetic averages of the corresponding (/, $o)~ 
elements of the respective generators. Denote EJ C the matrices obtained from Ef c by substituting the s-dependent 
terms with their averages: 



E fc = a/cEl s )( s + 2 l+^EI s )( s l+^EI s + 2 )(' 



a,f c h 2 



d fc I - 



bf c h~ 2 . 



(22) 
(23) 



The map given by the new Kraus generators E^ c is not trace preserving, so we are seeking for a set of generators 
Ef c which satisfy the trace preserving condition ^2f c El c Ef c = I. To this end we need to scale the circulant 

M, so the new trace preserving circulant Kraus generators are defined Ef c = 
M-V 2 . 



matrices Ef c by Y, f c E Jc E Jc 



(E fr E^E 
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'. + «2 h + k% h 



(24) 



where 

«i = E(|a/c| 2 + |d/c| 2 + |6/ c | 2 ) 

= j(>E^ <<&o|6j, +1 46i6 i+ i|* > + ^E;,y <*o|&S?&?|*o) + ^E^ <*o|&54 +1 &i+i&,|*o>) (25) 

«a = |(^E jy <*o|6^5, + i^-il*o> + fE w - <*o|65?6A+i|*o>). 

The matrix M is a circulant band matrix which can be diagonalized via discrete Fourier transform F (see Appendix 
B): 



fHif = V (/si + K 2 ^ 2 ' + k^- 2 ') |z)(z| 



i 



The matrix F'MF is a diagonal matrix with all diagonal values nonzero, so its inverse exists. The normalization 
operator A can now be defined as A = F{FHlF)- l l 2 F\ where (F^MF)' 1 / 2 = J2i ( K i + K 2^ 2 ' + n*uj- 2l y 1/2 \l){l\. 
A simple check verifies that A 2 M = / and A = AI~ X / 2 as desired. 

Unfortunately, the normalized generators Ef c — AEf c contain all even powers of the matrices h for general values 
of the Chern-Simons parameter k and do not admit a helpful form. However, for the special values k — 2, 4 there is a 
significant simplification. Notice that each of the terms in expressions for K\, H2 hi Eq. (25) is an averaged Kauffman 
bracket. For large N most of these brackets will correspond to disjoint links, i.e. the links formed by two step walks 
involving the strands located at j and f do not entangle. As shown in Table I, this always occurs if \f — j\ > 3. If 
we approximate the average by the value of the bracket for these disjoint links then 



K i - 32 



^(6cos(^)+4cos(^) + 2cos(^) + 5)sec 4 (^) + 0(l/iV) 

(26) 
«a = §(cos(^)-co S (^) + l) S ec 2 (^))+0(l/JV). 

where the error is of order 1/N. Thus at the special values k — 2,4 we have Ki ~ and the normalization operator 

— 1/2 

becomes a scalar multiple of the identity: A = n x ' I. Henceforth we focus on the case k — 2 since the behaviour of 
k = 4 is quite similar. The Kraus generators are given by 

Efc = n~ 1/2 (a fc h 2 + d fc T+b fc h- 2 ). (27) 

The elementary circular matrix h can now be diagonalized as F^hF = g~ = X) n ez w ™p)(^ where u> — e 2m ' N . This 
allows to write the Kraus generators in a diagonal form 

E fc = K- 1/2 F(a fc g 2 + d fc l + b fc g- 2 )F^ (28) 

= Y, X Mk)P fk (29) 

feeZjv 

with Xfc(k) = k[ (af c Lu 2k + df c + bf c uj~ 2k ) and Pf k = F\k)(k\F^. The action of the superoperator £ on a spatial 
density inatrix ps is then given by 

i( Ps )^J2 E ^P^l = J2 E *Mk)\%(l)P fh paP}, (30) 

f,c f,c kA^LltN 

and for t time steps 

£ t (p S )= E IlY<(. x f^ x U( l ))Pf k ps p l (si) 

fc.ZSZjv i=l fi,Ci 

Let's use this compact form to calculate the diagonal probability distribution p(s,t) at time t: p(s,t) — 
(s\£ t (ps(0))\s). Let's assume that the walker is initialized in the position eigenstate |sq), that the size of the 



periodic lattice is N, and that during the 2-step walk the coin is always reinitialized to state c = 0). We find for the 
index k — 2: 

Pising(s,t) = (s|£*(|s )(s |)|s) 

- 1 V uissoKr-Vf 2(r-l) , -2(r-l)\t 

- _J_V* MV a ,(8-«o)(r-0 w 2m(r-0 a ,-2(t-m)(r-0 
/(\P . .rfs-sn+4rn-2tl V^ -l(s-s a +4m-2t) V oz / 



1 V* fMV , ,r(s-s +4m-2t) y^ 

2' AT 2 Z^m=0 Vm/ Z^rSZjv W Z^ 

2^5" E m =0 ( m )(^^- Sn +4m-2t,o) 2 
1 ( * \ 



lei 



This is the binomial distribution where the range of the sites is s G [— 2£, 24] and the probabilities are nonzero only for 
s = so + 4n, n € Z, i.e. the y 2 model with fc = 2. Ising anyon walkers therefore have the same probability distribution 
as the classical random walk where every step moves two units to the right or left and the variance, scaled so that 
each two steps move takes place over two time intervals, is a 2 s in „(£) = t. 

IV. CONCLUSIONS 

We have considered a lossy anyonic quantum walk protocol where the entanglement between the spatial modes 
of the walker and its environment, the coin and fusion space, is lost on every second step. We calculated the time 
evolution of the exact probability distributions for various values of the parameter k in Chern-Simons theory. The case 
k = 1 corresponds to the Abelian anyonic quantum walk and it was found to be equal to the standard V 2 quantum 
walk with trivial exchange statistics, the variance of which has a leading term proportional to t 2 . Cases k > 2 are 
non- Abelian anyon models, for which the variance grows linearly as a function of time for all tested values of k in the 
time scales of up to 100 iterations of the superoperator. By approximating the Kraus generators of the superoperator 
by circulant matrices, the generators can be diagonalized for k — 2, 4. This allows for a compact expression for the 
probability distribution after arbitrary number of iterations of the superoperator. The expression for k = 2 is the 
binomial distribution which is equal to the classical random walk distribution. Thus the V 2 anyonic quantum walk 
with Ising anyons has the same behaviour as the fully coherent walk, which was also shown to have a linear variance 
with coefficient 1 [ ]. 

In the fully coherent walk, the slowdown of the walker propagation can be explained by decoherence: the fusion 
Hilbert space of non-Abelian anyons acts like an environment to the walker+coin system, degrading the quantum 
correlations between the spatial and coin modes which are the origin of the quantum speedup. It was conjectured in 
Ref. [ ] that the slowdown happens when the parameter k is much smaller than the number of time steps t, k <C t, 
ie. for any finite k, the walk behaves classically in the large time limit. In the V 2 model presented here the results 
support that conjecture. There are two kinds of decoherence mechanisms at work, one due to the fusion space of the 
anyons and the other because the entanglement between the spatial modes and the coin is lost. The loss of quantum 
correlations on every second step does not change the qualitative behaviour when the braiding statistics is Abelian, 
but when the additional effect of the fusion space comes into play for non-Abelian anyons, the spreading velocity 
changes to diffusive for all values of k for long enough times (in fact after fewer than 100 steps). 

The V 2 model can be easily generalized to V n models where the walk evolves coherently for n time steps instead of 
two. Evaluation of higher number of time steps becomes increasingly hard, because the number of different Kauffman 
brackets that need to be computed increases as the walker is allowed to do more steps between the tracing operation. 
One might expect that for large n, the walker propagates diffusively for n steps for small k and ballistically for large 
k, before the tracing is carried out. However, as we have shown, the variance of the anyonic V 2 walk is linear even 
when k 3> 2, so we expect that in the long time limit the variance is linear for all V n models, although the walker 
might spread ballistically in the initial stage of the walk. 

It is known that spatial randomness in the coin operator can lead to localization of the walker wave packet [13], ie. 
the probability to stay in the initial position becomes very high and the probability falls off exponentially away from the 
initial position. In the anyonic quantum walk setup, spatial randomness can also occur if the occupation numbers of 
the background anyons fluctuate randomly. These fluctuations might occur when thermal excitations create particle- 
antiparticle pairs from vacuum, for example. It is interesting to note that if the occupation numbers of the islands 
are not uniform, the Abelian anyons can pick up different phases during forward and backward time evolution, and 
the effect of the braiding statistics becomes nontrivial. We have calculated evolution of the V 2 anyonic walk with 
random fillings of Abelian anyons on the islands, but observed diffusive spreading (no localization). Localization is a 
quantum phenomenon that requires interference of the probability amplitudes to occur. Interestingly our results show 
that while the memory loss in the V 2 model preserves enough coherence in the Abelian walk to provide for quadratic 
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speed up without disorder, it does not provide enough to give localization with disorder. Rather the competition of 
the two effects yields classical behaviour. 
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Appendix A: Braid group matrix elements 



The construction of the superoperator £ requires calculating the eight expectation values in Eq. (8) for all values 
of s and s' . The expectation values are related to the Kauffman bracket via Eq. (5). Kauffman bracket polynomial 
(l/)(A) in the variable A [10] is a link invariant for framed, unoriented links. The defining properties of the Kauffman 
bracket are given by the uncrossing relation and two relations for removing loops from the bracket: 



= A 



A~ 



(Al) 



(LUQ) = -(A 2 + A- 2 )(L) 



(A2) 



(0) = L 



(A3) 



The value of the Kauffman bracket for a specific braid word can be calculated by forming the link diagram that 
corresponds to the braid presentation and Markov closure of the word. The unique link diagrams for one braid word 
are drawn for illustration in Fig. 4. The value of the bracket polynomial is then obtained by applying the relations 
(Al) and (A2) to the link diagram, removing all crossings and loops in the diagram except one. The bracket of a 
single loop is trivial, so the remaining coefficient is the value of the bracket polynomial. 





FIG. 4: Two links corresponding to the expectation value ( < &o|frs/&s/ +1 bs-i|<l?o)- a) Case s' — s — 1. The forward braids &s_! 
and the backward braids b s ,b s , +1 act on separate sets of strands, therefore the links corresponding to the forward and backward 
braids are disjoint. The links corresponding to s' — s < —4, s' — S > 1 are all disjoint and the values of the Kauffman brackets 
are the same, b) Case s' — s — — 1. The forward and backward braids now form a joint link. 

For SU(2)fc anyons, the value of the parameter is A = ie~ 47I 7 2 ( fe + 2 ) ; and the quantum dimension satisfies d = 
— (A 2 + A~ 2 ) — 2 cos( jtztj). For s' — s large, the forward and backward braid words never touch the same strands, the 
links corresponding to forward and back evolution are disjoint, and the value of the Kauffman bracket polynomial is 
equal for all s and s' . Thus, the calculation of a general element ^<&o |i?(s, s')|$o) involves only the calculation of the 
disjoint element and a few cases where s' — s is small. The unique values of five of the expectation values are given 
in Table I. For the remaining elements, we have 



(*o|#_i&2-i|*o> = (*o|&if&2|*o) 



t2j.2| 
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sf — s 


(*0 b\, 


_ 1 6t,_ 2 6 s _ 2 6 s . 


-1 $o) 


($o bl,bl +1 b s+1 b s $„) 


-2 


d~ 4 


d" 4 


-1 


d~ 2 


d" 2 





1 


1 


1 


d~ 2 


d" 2 


2 


d- 4 


d~ 4 


< -3, > 3 


d~ 4 


d- 4 



s/ — s 


(®o\bl*b s -2b s -i\G> ) 


<*o|&J,&J, +1 tf-i|*o) 


-3 


-A 6 (A 4 + A~ 4 )/d 3 


-A- 6 (A 4 + A- 4 )/d 3 


-2 


d- 2 


d- 2 


-1 


d- 2 


d- 2 





- J 4 6 (4 4 + J 4- 4 )/d 3 


-A- 6 (A 4 + A- 4 )/d 3 


< -4, > 1 


-A 6 (A 4 + A~ 4 )/d 3 


-A- 6 (A 4 + A- 4 )/d a 



sf — s 


<$o|&!?& 2 |$o) 


-1 


(A 4 + A" 4 ) 2 /d 2 





1 


1 


(A 4 + A~ 4 ) 2 jd 2 


< -2, > 2 


(A 4 + A- 4 ) 2 /d 2 



TABLE I: The values of the braid group elements involved in the evolution CP map. 
and the conjugate transpose elements 

(*0|&J,_1&I'-2*2|*0) = ((*0|&!V- a &y-l|*0»* 

and 

<#o|&I?_ 1 W.|*o) = «*o|6l6l+i^_i|*o»*. 



Appendix B: Circulant matrices and CP map generators 



A matrix C = (cy) of order n is called a circulant matrix [14] if Cj 



2 i-j( mod «)• The entries of the first column 
a = (ao, ai, ...,a n _i) determine the entire circulant matrix which we denote C n — circ(a) = circ(ao,ai, ...,a n _i), 

and in matrix form it reads 

Q-0 ffln-l On-2 • ' ' &1 

ai ao a n -i ■ • • 02 
C„ = a 2 a • • • a 3 (Bl) 

\ a n _i a„_2 a„_3 ■ • • ao / 



Alternatively circulant matrices are written in terms of elementary circulant matrix h = X^me-fo 1 n-i> \ m /\ m + 
1 = circ(0, 1, ..., 0), as C n = circ(ao, a%, ..., a n —i) = p(h) where p is the polynomial p(z) := a + <X\Z + ■ ■ ■ a n —\Z 
Matrix h is diagonalized by the finite Fourier transform unitary matrix F, with elements F a b 



n-l 



1 , ,ab 



e i2ir/n^ ag pjfop — g^ w here g = diag(l,w, ...,u> n x ). Any circulant matrix is then canonically decomposed as 
C n — circ(ao, ai, ...,a n _i) = F'diag{p{l),p(oS), ...,p(uj n ~ 1 ))F. A banded circulant matrix is a circulant matrix 
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for which only a connected subset in the sequence a — (aj)^Z is non zero. After introducing the approxi- 
mations issued in Eqs. (16)-(18) the generators Et c of QW's CP map become banded circulant matrices i.e. 
Ef c = circ(a = d fc ,ai — 0, a 2 — b fc , ...,a n _ 2 = a/ c , a„_i = 0). 
Circulant matrix h = X) s gz I s ) ( s + -M ' an( i its inverse h~ x 



Z>+i> 



= 5>-l)(n| = 






1 \ 


1 




1 




1 




V 


i o) 


/o 1 


\ 


1 







1 




'•• 1 


\1 


oy 



(B2) 



(B3) 



generate the abelian group {h a } a( zz N — I*n, share the property hf = h 1 = J2sei I s + ■"•) ( s l > where ^ 2 = 1- 

An important tool used in Sec. Ill to approximate the V 2 quantum walk model is the optimal circulant of a matrix. 
For an arbitrary square matrix D = (djk),i,k = 0, 1, ...,N — 1, we choose a suitable circulant matrix C. For the 
construction of C it has been recommended [I r >], (after preliminary transformations such as changing the order of the 
rows and multiplying them with suitable constants), to apply the following: by imposing periodicity i.e. dj+N,k = dj t k, 
form the arithmetic averages of the elements along the diagonals i.e. a,j — -h X^fcJo ^j+k,ki an d construct in this way 
the circulant C = circ(ao,ai, ...mazv-i)- 

On the other hand an optimal circulant approximation of some Toeplitz matrix T (see e.g. Strang's suggestion in 
[16]), have been constructed by determing the nearest circulant with respect to Frobenius matrix norm for matrix 
T, i.e. by solving the optimization problem mine -.circulant \\T — C\\ F [ ]. The resulting circulant coincides with the 
one obtained by the method of arithmetic averages [15], applied to T. In the quantum walk context the method of 
arithmetic averages is been applied to the Kraus generators Ef c that are approximated by the optimal circulant Ef c . 
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